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Abstract: For almost a century, the cosmological constant has been a 
mysterious object, in relation to both its origin and its very small value. By 
using a Bose-Einstein condensate analogue model for gravitational dynam- 
ics, we address here the cosmological constant issue from an analogue gravity 
standpoint. Starting from the fundamental equations describing a system 
of condensed bosons, we highlight the presence of a vacuum source term for 
the analogue gravitational field, playing the role of a cosmological constant. 
In this simple system it is possible to compute from scratch the value of this 
constant, to compare it with other characteristic energy scales and hence 
address the problem of its magnitude within this framework, suggesting a 
different path for the solution of this longstanding puzzle. We find that, 
even though this constant term is related with quantum vacuum effects, it 
is not immediately related to the ground state energy of the condensate. On 
the gravity side this result suggests that the interpretation and computation 
of the cosmological term as a form of renormalized vacuum energy might be 
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misleading, its origin being related to the mechanism that instead produces 
spacetime from its pregeometric progenitor, shedding a different light on 
the subject and at the same time suggesting a potentially relevant role of 
analogue models in the understanding of quantum gravity. 

0.1 Introduction 

The cosmological constant pQ has been one of the most mysterious and 
fascinating objects for both cosmologist and theoretical physicists since its 
introduction almost one century ago [2] . Once called by Einstein his greatest 
blunder, it seems nowadays the driving force behind the current accelerated 
expansion of the universe. At the theoretical level, the explanation of its ori- 
gin is considered one of the most fundamental issues for our comprehension 
of general relativity (GR) and quantum field theory (QFT). 

Since this constant appears in Einstein's equations as a source term 
present even in the absence of matter and with all the symmetries of the 
vacuum (i.e. a stress-energy tensor of the form T^ v oc g^ u ), it is usually 
interpreted as a vacuum energy, essentially related to the zero point fluctua- 
tions of quantum fields. This reasonable point of view has originated what is 
often called the "worst prediction of physics" . Indeed, its theoretical value, 
which is naively obtained by integrating the zero-point energies of modes 
of quantum fields below Planck energy (but can be computed also by more 
sophisticated renormalization arguments), is about 120 orders of magnitude 
larger than the measured value. We can summarize the situation by saying 
that, given the absence of custodial symmetries protecting the cosmological 
term from large renormalization effects, the only option we have to explain 
observations is fine tuning O 0] . 

Despite the large number of attempts (most notably supersymmetry (5], 
which, however, must be broken at low energy) to improve this estimate, 
this problem is still waiting for a complete and satisfactory solution. 

This huge discrepancy is plausibly due to the use of effective field the- 
ory (EFT) calculations for a quantity which can be computed only within 
a full quantum gravity (QG) theory (see, however, [6] for a proposal in the 
semiclassical gravity limit). Unfortunately, to date, we do not have any 
conclusive theory at our disposal, allowing a complete calculation of the 
gravitational effective action from basic principles, and solving the natural- 
ness problem associated to the cosmological constant. 

Therefore, in order to give support to this idea, some other arguments 
have to be presented. Analogue models for gravity [7], in this respect, give us 
the possibility to study the emergence of given EFT with certain geometrical 
content from microscopic constitutents (typically interacting atoms, fluids, 
etc.), keeping the process under control at every stage of the transition, 
and hence allowing us to show at which point of the calculation the EFT 
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intuition of the cosmological constant as a vacuum energy term fails. 

Already in [HJ [9j [10] it was shown that a naive computation of the ground 
state energy using the EFT (the analogue that one would do to compute the 
cosmological constant) would produce a wrong result. The unique way to 
compute the correct value seems to start from the full microscopic theory, 
working out the macroscopic quantities from it. 

Given the deep difference in the structure of the equations of fluid dy- 
namics and those of GR (and other gravitational theories) it is not possible 
to have an accurate analogy at the dynamical level: indeed, this is forbid- 
den by the absence of diffeomorphism invariance. However, in [11] it has 
been shown for the first time that the evolution of the acoustic metric in a 
Bose-Einstein condensate (BEC) is described by a Poisson equation for a 
nonrelativistic gravitational field, thus realizing a (partial) dynamical anal- 
ogy with Newtonian gravity. Noticeably, this equation is endowed with a 
source term which is present even in the absence of real phonons and can 
be naturally identified cosmological constant. 

In this chapter we will consider such analogue model for gravity and 
directly show that the cosmological constant term cannot be computed 
through the standard EFT approach [T2], confirming the conjecture of OS]. 
Somehow unexpectedly, we find that also the total ground state energy of 
the condensate does not give the correct result: indeed, the cosmological 
constant is comparable with that fraction of the ground state energy corre- 
sponding to the quantum depletion of the condensate, i.e. to the fraction of 
atoms inevitably occupying excited states of the single particle Hamiltonian. 

This result is twofold. First, it gives an explicit calculation showing 
qualitatively and quantitatively where the EFT intuition might fail to grasp 
the nature of the vacuum term. Second, it shows that the subject of analogue 
models is not just a mere curiosity for condensed matter physics, but it is 
also of great interest for research in quantum gravity, for their ability in 
providing guidance and patterns to be used to address what is the most 
urgent problem there, that is the recovery of the continuum semiclassical 
limit. 

0.2 The cosmological constant problem in EFT 

The cosmological constant A enters Einstein's equation 



as a term multiplying the metric tensor g^ v . This means that we can inter- 
pret A as a vacuum energy density 
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and hence see the cosmological term as a vacuum stress energy term 

T^ = -£g, v . (3) 

The other way around, this argument implies that every vacuum energy 
is in the form of a cosmological constant term and must be thus identified 
with A. 

The problem arises when computing the renormalization of this term 
using quantum field theory techniques. Summing up all the zero-point en- 
ergies of all normal modes of some quantum field of mass m up to a cutoff 
energy /i, one obtains [T3] 



r lhc 4tt k 2 dk 1 f~ m 2 c 2 ^ 
£ = J 1^2W fc+ ^"i^W () 

If we compute this term with a Planck-scale cutoff, we obtain [1] an estimate 
for the vacuum energy of the order of 

£ P » 10 110 erg/cm 3 , (5) 

while if we lower the cutoff to a much smaller energy scale, like the QCD 
scale, we expect 

£qcd~ 10 36 erg/cm 3 . (6) 
The cosmological observations give 

Sobs- 10" 10 erg/cm 3 , (7) 

a value hardly explained in terms of obvious particle physics scales. 

In a proper semiclassical renormalization procedure of the whole theory 
(gravity plus quantum fields) a suitable bare cosmological constant Ab has 
to be introduced to tune the renormalized A to match its observed value p3] . 
However, since the order of magnitudes of Aeft and the observed one are 
so different, such a procedure would require an extremely fine tuning of Ab- 

Furthermore, another obscure issue should be clarified. In quantum 
field theory in Minkowski spacetime, the zero-point energy has no particular 
meaning, since energy does not gravitate. Actually, it can be safely removed 
by normal ordering. On the contrary, when gravity is turned on, this is no 
longer possible, since energy directly enters Einstein's equations: besides 
the very difficulty of defining a normal ordering prescription in curved dy- 
namical spacetimes, different zero-point energies yield different stress-energy 
tensors T^v and, therefore, different solutions of Einstein's equations. 

What is then the correct vacuum energy that must be put into Einstein's 
equations? 

From such considerations, it seems that there is no way out to this 
problem. The only possibility to compute A would be to directly know 
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the underlying QG theory from which both GR and EFT emerge as a low 
energy limit. Under this hypothesis there is no hope of computing A from 
scratch working only in a semiclassical gravity framework. In a certain sense 
A is an emergent quantity that parametrizes in one number certain details 
of the microscopic structure and dynamics of spacetime and influences the 
macroscopic equations of GR and EFT, and, as such, cannot be justified 
avoiding a discussion of the microphysics. 

As an example of this emergence mechanism, let us mention the propa- 
gation of phonons in BECs. In that case, the equation of motion of phonons 
can be written using only microscopic quantities (speed of sound c, velocity 
v, and density of the fluid p) plus the healing length 



which, through the combination of microscopic quantities (m is the mass of 
bosons), defines the scale at which the dispersion relation 



is no more in the relativistic regime oj 2 = c 2 k 2 . If an observer could only 
make measures of phonons, he/she would measure the healing length, but 
he/she would not be able to compute it from first principles, just because 
he/she would not know the whole theory of BECs but only the phonon EFT. 

To shed some light on this problem, it is interesting to have a toy model of 
gravity where one can compute from first principles both the vacuum energy 
and the cosmological constant. In this way one can compare them and check 
if they are the same quantity or they are unrelated instead. In order to do 
this, one needs a model where it is possible to derive not only the dynamics 
of a field leaving in an effective geometry described by some metric g^, as 
usual in analogue models [7] , but also the equation governing the dynamics of 
the metric itself. That is, analogue Einstein's equations are needed, so that 
the analogue cosmological constant can be directly read from them. Such a 
model was studied in [11], using a BEC with C r (l)-symmetry breaking. The 
computation of the analogue cosmological constant is performed in [T2] and 
reported in Sec. 10.41 together with a revision of the model of [11] . 

0.3 Volovik's proposal 

The idea of using analogue models to understand the origin of the cosmolog- 
ical constant is not original of [12] but was firstly developed by Volovik [8j [9] 
for a quantum liquid and presented in [10] also for a Bose gas. However, 
his approach is completely different from that of [12]. It is worth briefly 
reviewing it, to make a comparison with ours. 



n 



(8) 





(9) 
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Even if Volovik did not have an analogue model with an equation de- 
scribing the dynamics of the geometry as in [11] , he nevertheless used a nice 
argument to determine the analogue cosmological constant. His argument is 
based on the identification of the proper thermodynamical potential for the 
particular considered problem. In this case, we are interested in the emer- 
gence of an analogue QFT in condensed matter. The many-body system of 
identical atoms constituting the quantum liquid is described by the grand 
canonical Hamiltonian 

H = H- fiN, (10) 

where H is the second-quantized Hamiltonian, /i is the chemical potential, 
and TV" the particle number operator. 

The correct vacuum energy density for the QFT emerging in the many- 
body system (corresponding to the analogue cosmological constant) is there- 
fore the expectation value of Eq. ()10p on a state |0) with no phonons, in the 
thermodynamic limit in which both the volume V and the particle number 
N goes to infinity 

£ vac = ^(0|#-/xiV|0>. (11) 

Using the Gibbs-Duhem relation of thermodynamics [15], stating that at 
thermodynamic equilibrium 

E — TS — fiN = -pV, (12) 

where E and N are the expectation value of H and N, respectively, T is the 
temperature of the system, S its entropy, and p its pressure, one obtains 

£vac = ~P, (13) 

because T = in the phonon ground state |0). 

Equation (|13p is the key result of this analysis, since it exactly represents 
the correct equation of state for the cosmological constant. To summarize, 
if the vacuum energy is the density of the expectation value of the grand 
canonical Hamiltonian on the zero-temperature state (no excitations), then 
p = —£ vac by thermodynamic relations. The second interesting feature is 
that, if % = (H) of Eq. ([10")) were really the quantity that gravitates in 
place of the energy E = {H), the freedom in the choice of the zero-point 
energy would not affect the gravitating quantity T~L. If the energy of each 
atom were shifted by a factor of a, the Hamiltonian would be shifted by 
aN . However, also the chemical potential would have to be shifted of a, 
such that the grand canonical Hamiltonian, being the difference of H and 
[iN, would not change under this transformation. In so doing, there would 
be a definite gravitating quantity, allowing at the same time for the freedom 
in the choice of the zero-point energy. 

Unfortunately, this argument does not use any dynamical equations for 
the spacetime geometry. This is, in our opinion, the weakness of the above 
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treatment. Indeed, the cosmological constant is the quantity that gravi- 
tates in the absence of matter, encoding the microscopic properties of the 
spacetime structure. It may be a vacuum energy of some field, but this 
interpretation may be also wrong. As discussed in the previous section, the 
only way to obtain a reliable result for the cosmological constant would be to 
derive a dynamical equation for the metric, from which A might be directly 
read off. 



0.4 A lesson from BECs 

As we have argued above, the only safe way to compute the cosmological 
constant would be to know the microscopic structure of the spacetime and to 
derive Einstein's equations from it. In this section we apply this procedure 
to a particular analogue system, a BEC with U{1) breaking [11], for which 
an equation describing the analogue gravitational dynamics exists. 

0.4.1 Settings 

The model used in [11] is a modified BEC system including a soft breaking 
of the U(l) symmetry associated with the conservation of particle number. 
This unusual choice is a simple trick to give mass to quasiparticles that are 
otherwise massless by Goldstone's theorem. In second quantization, such a 
system is described by a canonical field satisfying 

[%x),#t( t y)] = ( 5(3)( x _ x >) 5 (14) 

whose dynamics is generated by the grand canonical Hamiltonian 

n = H-fiN, (15) 



where 

H = Jd 3 : 



T? A 

* t V7 \T/ _l T/vf/t\f/ _i_ ^\T/t>f/t\f/\f/ I \f/\f/ _i_ if/tif/t 



V X ^ T V x tf + VW^> + — \I/T\|>T\|/\|/ _ _ ^ + \J/T\J/ 

2m 2 2 



(16) 

is the Hamiltonian and 

N = [ d 3 x¥^ (17) 



is the standard particle number operator for In the Hamiltonian, g is the 
coupling constant of the two-body interaction, while A, having dimensions 
of energy, represent a U(l) breaking term, associated to a violation of the 
conservation of the number operator. For further details on this model and 
on possible physical realizations, see [11]. See also [16] for a generalization 
to condensates with many components. 
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We describe the formation of a BEC at low temperature through the 
complex function f° r the condensate and the operator <f) for the pertur- 
bations on top of it, defined by 



* = * (1 + 



(18) 



The canonical commutation relation is directly obtained from the commu- 
tation rules of the boson field ^ [Eq. (|14p]: 



[0(t,x),^(t, x ')] = -^T <^ (3) (x-x'). 

pom 



(19) 



Using Eq. (JTHJ) , it is convenient to expand the grand canonical Hamiltonian 
H of Eq. (fT5]) up to second order in <j> 

H^H + Hi+H 2 , (20) 



where 
■H = J d 3 x 
Hi = I d 3 : 

Ho 



-^Vl + V-»+Zpo)y -*(y 2 o + n 2 ) , (21) 



2m 



2m 



Vi + V - {i + gpo ] *o - A*f f 



+ h.c, (22) 



/ 



d 3 xp 



ft 2 V^r 



p 2m f f 



+ V - n + 2gp 



+f (0 t2 + 2 ) + ^(^ 2 + ^ t2 )}- (23) 



For a stationary condensate, dt^fo = and the grand canonical Hamil- 
tonian ([22]) generates a modified Gross-Pitaevskii equation 



h 2 

2m Wq 



*n = 0. 



(24) 



Moreover, to compute the analogue cosmological constant, it is enough to 
consider homogeneous backgrounds. Thus, one can assume that V = and 
the condensate is at rest, such that has a constant phase, that one can 
put to (\&q = = y/po)- With these assumptions, Eq. (f24"|) simplifies to 



p = gpo - A. 



(25) 



Under the same assumptions, the equation of motion of the quasiparticles, 
generated by the second order Hamiltonian (|23p . reads 



ihd t 



'^i + gp + \ 
2m 



(26) 
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To solve this equation it is convenient to define at first a two-component 
field pi] 



w (27> 

Then Eq. (|26p can be written in a compact form 

ihd t W = BW, (28) 
B = (T + gp + X)a 3 + i{gpo-X)a 2 , (29) 

where T is the kinetic energy operator 
and <Tj are the Pauli matrices 

-(?;)■ -(i-0- < 3i > 

Since the field W is invariant under the conjugation operation defined by 

S = ffiS*, (32) 

the structure of must be 

W = J d 3 k(WA + W k a{), (33) 

where W k is a doublet of C-functions. 

Using Vj = — V x , T T = T and the properties of Pauli matrices, one 
verifies that the scalar product 

{W!\W 2 ) = Jd 3 x Po (x)Wl T (t,x)a 3 W 2 (t,x) (34) 

is conserved under time evolution when Wj are solution of (|28|) . since 

B* T a 3 = a 3 B. (35) 

Imposing the following normalization for the modes 

<W k |W k /) = -{W^Wv) =^ 3 )(k-k'), (36) 

(W k |W- k /) = 0, (37) 

one gets 

[6k, al] = [{W k \W},-{Wv\W)] = (W k |W k .) = ^ 3 )(k-k / ), (38) 
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which shows that and are in fact destruction and creation opera- 
tors. Moreover, all the other scalar products and all the other commutators 
vanish. 

By homogeneity and stationarity, the doublets must have the following 
form 



-iurf+ik-x 



(39) 



where y pq(2tt) 3 is a convenient normalization factor and and are 
constant. The field (p can now be expanded as 

d 3 k 



-iajt+ik-x ~ 



a k + v k e 



* +iwt— ik-x ~ t 



(40) 



and the Bogoliubov coefficients itk and i>k obey the standard normalization 

1, (41) 



i 1 9 i i2 
|«k| - Pk| 



which follows by inserting the above expansion of <j) m Eq. (|36|) . In term of 
these Fourier components, Eq. (|26p reads 

+ gpo + A 



— 



2m 

h 2 k 2 

2m 



+ 5Po + A 



«k = (5Po - A) v k , 

- {gpo - A)« k . 



(42) 



i'k 



By imposing that this system has nontrivial solutions, one obtains the quar- 
tic dispersion relation 



h 2 oo 2 = 4Xgpo + 



4 1,4 



9Po + X h 2 k 2 1 



m Am 2 ' 

describing massive phonons with ultraviolet corrections, mass M 

M = -m, 



and speed of sound c s [TT 



9Po + A 
2 ffPo + A 



m 



(43) 



(44) 



(45) 



Finally, by using the mode normalization (|4ip . the system (|42p and the 
dispersion relation (j4"3|) . it is possible to compute «k and v^: 



l-D 2 k 



l-Dl 



where both u k and t>k are chosen to be real and 

Huj - (fr?k 2 /2m + gp + A) 



gpo - A 



(46) 



(47) 
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0.4.2 Analogue gravitational dynamics 

When the homogeneous condensate background is perturbed by small in- 
homogeneities, the Hamiltonian for the quasi-particles can be written as 

(seem 

i^quasip. « Mcl - + M$ g . (48) 

-f^quasip. is the non-relativistic Hamiltonian for particles of mass M [see 



Eq. (|44p ] in a gravitational potential 

».(x) = (9P0 + 3 2 ^ P0 + A) VW (49) 

and 8p(x) = [(/9o(x)/poo) — l]/2, where is the asymptotic density of the 
condensate. Moreover, the dynamics of the potential $ g is described by a 
Poisson-like equation 



x L 2 



$ g = 4vrG 7V Pp + C A , (50) 
which is the equation for a non-relativistic short-range field with length scale 



L= - = (51) 
Vievrpoa 3 



and gravitational constant 



r = 9(9 Po + 3A)(gp + A) 2 
N 4^ 2 mA 3 /2( OT )i/2 ' 

Despite the obvious difference between <£ g and the usual Newtonian grav- 
itational potential, we insist in calling it the Newtonian potential because 
it enters the acoustic metric exactly as the Newtonian potential enters the 
metric tensor in the Newtonian limit of GR. 

The source term in Eq. (|50p contains both the contribution of real 
phonons (playing the role of matter) 



(CWlO - {o\ftm) + ^Re f(CI#IC) - <0|#|0) 



(53) 



where |£) is some state of real phonons and |0) is the Fock vacuum of the 
quasiparticles (ak|0) = 0, Vk), as well as a cosmological constant like term 
(present even in the absence of phonons/matter) 



n Zgpo (gpo + 3 A) (gp + A) 
C A = ^ Re 



(O|0t0| O ) + I(O|#|O) 



(54) 



Note that the source term in the correct weak field approximation of Ein- 
stein's equations is AttGn(p + 3p/c 2 ). For standard nonrelativistic matter, 
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p/c 2 is usually negligible with respect to p. However, it cannot be neglected 
for the cosmological constant, since pa/c 2 = — p\. As a consequence the 
analogue cosmological constant is 

0.4.3 BEC ground state energy 

By using the formalism developed in Sec. 10.4,11 we compute the vacuum 
expectation value of H in the ground state |0). To this aim, it is convenient 
to use the expansion of % in powers of <j) given in Eq. (|20p . The energy 
density h$ of the condensate (density of Hq) in an homogeneous condensate 
is straightforwardly obtained from Eq. (|2ip . using the relation (|25p between 
the chemical potential /i and the couplings g and A 

h = - 9 -^. (56) 

The expectation value of jii vanishes because it contains only odd powers 
of the phonon field cj). Finally, the density hi of the expectation value of H.2 
can be computed by using the equation of motion of the perturbations ([26]) . 
such that the expression of %2 simplifies to 



K2 = j /d 3 x/ 



ftdtt - . (57) 



Inserting the field expansions ()40p in the above expression, one obtains 



(2k) 



(58) 



which implies that the density of (0|7^2|0) is 



/c\ 3 k 
J^^\ V]l \\ (59) 

The integral in Eq. ([59]) is computed by using Eqs. (f4U|) and (|47p. The 
divergence in d = 3 spatial dimensions is regularized by performing the cal- 
culation with d < 3 and then going to the limit d — )■ 3. This regularization is 
equivalent to the subtraction of higher order interaction terms (see also [15] ) 



h 2 = Mgfiyfcffi F h (±) , (60) 

15 V 71 " \9PoJ 

where a = 4irgm/h 2 is the scattering length, is plotted in Fig. Q] (dashed 
line) and -FX(0)=1. The total grand canonical energy density is therefore 
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Figure 1: F h [dashed line, Eq. (|fiT)|)]. F p [dotted line, Eq. fBT])]. [dot- 
dashed line, Eq. ([70])]. and F A [solid line, Eq. (fTTj) ]. 



h = ho+h2 = £ii r 



128 > =r r / A \ 

15y 7r V3/W 



(61) 



To compute the energy density e = h+fip, that is the density of (0|i?|0) = 
(0\T-L + pN\0)), we have to express at first the density of condensed atoms po 
in terms of the total number density p. To this aim, we expand the particle 
number operator N in powers of <j) 



N = N + N 1 + N 2 



(62) 



where, as in Eq. (|2Up . Nq, N%, and N 2 contain respectively no power of 
only first powers, and only second powers 



N 
N 2 



d 3 xp , 

d 3 x $ + h.c. 



The density of Nq is then simply 

Po = l^ol 2 , 

(0|iVi|0) vanishes, and the density of (0|iY 2 |0) is 



(2^)3 



\ v k\ 



8po 
3^ 



\9PoJ 



(63) 
(64) 
(65) 

(66) 
(67) 
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where F p satisfies F p (0) = 1 (see Fig. [I] dotted line). This is the number 
density of non-condensed atoms {depletion). Note that pga 3 is the so called 
dilution factor which has to be much smaller than 1 for the Hamiltonian (|16p 
to hold. 

Furthermore, when A = 0, inverting the expression for total particle 
density, p = po + P2, one obtains, up to the first order in \J pa? 

(68) 

which is the density of condensed atoms in terms of the total density p and 
the scattering length a [IS] . In this case, p = gpo, such that the energy 
density is 

9P 2 

e = h + up = — 

This is the well known Lee-Huang- Yang [18] formula for the ground state 
energy in a condensate at zero temperature. In general, when the U(l) 
breaking term is small, this term is expected to be the dominant contribution 
to the ground state energy of the condensate. 



Po 



1 + 



128 

I5v/? 1 



pa^ 



(69) 



0.4.4 What does the cosmological constant correspond to? 

We shall now compare the energy density and the grand canonical energy 
density found in the previous section with the effective cosmological constant 
A of Eq. ([55]) . Ca of Eq. (|54p is computed by using Eq. ([6"T]) and the following 
expectation value 

/d^/c 8 / A \ 
— — 3 u k v k = -= \/^3 F I ( 70 ) 

where F^{Q) = 1 (see Fig. [H dot-dashed line). We finally obtain 

A = _ 20m^ + 3A ) v ^ / M 
3y/Trh 2 X \9PoJ 

where F A = (2F p + 3F^)/5 (see Fig. [Q solid line). 

Let us now compare the value of A either with the ground-state grand 
canonical energy density h [Eq. ([6"T])]. which was suggested in [SJ [9] as the 
correct vacuum energy corresponding to the cosmological constant, or to the 
ground-state energy density e of Eq. (|69p . Evidently, A does not correspond 
to either of them: even when taking into account the correct behavior at 
small scales, the vacuum energy computed with the phonon EFT does not 
lead to the correct value of the cosmological constant appearing in Eq. (|50p . 
Noticeably, since A is proportional to \J poa 3 , it can even be arbitrarily 
smaller both than h and than e, if the condensate is very dilute. Further- 
more, A is proportional only to the subdominant second order correction of 
h or e, which is strictly related to the depletion [see Eq. (f67l) ]. 
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Furthermore, several scales show up in the emergent system, in addition 
to the naive Planck scale computed by combining the emergent constants 
Gn, c s and h: 

L ^\S K (^)" 3/4( "° o3) " i/4a ' <72> 

For instance, the Lorentz- violation scale (i.e., the healing length of the con- 
densate) 

Lw = £ oc (p « 3 )~ 1/2 « (73) 

differs from Lp , suggesting that the breaking of the Lorentz symmetry might 
be expected at a much longer scale than the Planck length (much smaller en- 
ergy than the Planck energy), since the ratio Lvj/Lp oc (po<2 3 ) _1//4 increases 
with the diluteness of the condensate. 

To conclude, it is instructive to compare the energy density correspond- 
ing to A 

Ac 4 

£ * = «fc < 74 > 

to the naive Planck energy density 

fp = i- < 75 > 

The former is much smaller than the value computed from zero-point-energy 
calculations with a cut off at the Planck scale. Indeed, the ratio between 
these two quantities 

^oc Po a 3 (— V 5/2 (76) 
£p \9PoJ 

is again controlled by the diluteness parameter poa 3 . 



0.4.5 The spinor BEC case 

From the previous analysis one can easily realize that Llv scales with po fl3 
exactly as the range of the gravitational force [see Eq. (|5ip], signaling that 
this model is too simple to correctly grasp all the desired features. How- 
ever, in more complicated systems [16], this pathology can be cured, in the 
presence of suitable symmetries, leading to long range potentials. 

In the case of a BEC system containing several different bosonic species, 
the so called spinor BEC (see in these proceedings for a discussion of 
their relevance as analogue models and for references) it is possible to re- 
peat the analysis that we have described here, with minimal variations in 
the method. There are significant differences, however. In general, without 
tuning or symmetries (internal symmetries among the various components, 
essentially) the geometrical structure describing the propagation of phonons 
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is not a Lorentzian metric, at low energy, but rather something like a Finsler 
structure [20]. Additionally, the coupling between phonons and the conden- 
sate wavefunctions becomes extremely nontrivial, with the analogue New- 
tonian potential obeying some complicated equation that cannot be cast in 
the form of a Poisson equation, even including a Yukawa mass term. See 
[TC] for the complete discussion. 

However, if one imposes a symmetry in the system, that is, if there is 
an underlying symmetry under the permutation of the species, the situation 
improves so much that indeed a realistic analogue can be obtained. Con- 
sider a system with N components (whose nature will be neglected, for the 
present reasoning), labelled by roman letters A, B = 1, N. A Hamiltonian 
obeying the requirement just mentioned would be of the form: 

H= [d 3 xjr* A V(x) + f |^| 2 + { £ \*bA * A 

J A=\ \ B^A J 

(77) 

where /x, A, g are playing the same role of the corresponding quantities de- 
fined in ()16p . In this case, due to the richer structure available, one can 
recover a low energy notion of Lorentz invariance, and a distinguished long 
range analogue Newtonian potential. 

In particular, one can describe the deviations from homogeneity in each 
component of the condensate by expanding the wave function as 

[*oU = V[Po]a + a A {x) + tf A {x) , (78) 

where a(x),/3(x) are real functions. Within this ansatz it is then possible 
to show (see again [16]) that a long range Newtonian potential of the form 

N 

^ n {x)oc^2Pa(x). (79) 
A=l 

can exist and that this potential is indeed coupled to several (massive) 
phonons in a universal way, hence manifesting the emergence of an ana- 
logue of the equivalence principle. 

Therefore, the fact that in the single BEC there seems to be a basic 
flaw given by the short range potential should be seen as an artifact of 
the simplicity of the model, and not a basic obstruction. The above men- 
tioned generalized model can remove this unpleasant feature, allowing the 
construction of a more realistic dynamical analogue. 

Most importantly, the very same analysis that has been reported here 
can be done for the multi-component case so to study the nature of the 
vacuum contribution to the source term of the analogue Poisson equation. 
Even without reporting a full calculation one can easily see that the bulk of 
the result will be of the same nature: the cosmological term will be related 
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to the depletion factor, and hence non necessarily Planckian (in the sense 
discussed above). 

0.5 Summary and Conclusions 

Let us summarize the results. We have investigated the problem of calculat- 
ing the cosmological constant through an analogy with condensed matter. 
As a model we adopted the BEC with [/(l)-symmetry breaking proposed 
in |11] , given that in this system it is possible to define the analogue of the 
gravitational field satisfying a modified Poisson-like equation. The source 
term appearing in this equation is made of two pieces: the first one is re- 
lated to real phonons, corresponding to matter fields, in the analogy, and 
that generates the gravitational field associated to the presence of clumps 
of matter. The second one represents a vacuum contribution, that precisely 
matches a cosmological constant in the Newtonian regime. 

Of course this is just a toy model for gravity, which nevertheless provides 
a clear description of the procedure that should be followed to correctly com- 
pute the cosmological constant. We showed that the analogue cosmological 
constant cannot be computed as the total zero-point energy of the condensed 
matter system, even when taking into account the natural cut-off coming 
from the knowledge of the correct microphysics, as suggested in [81 [9]. In 
fact, the value of A is related only to a part of the zero-point energy, namely 
a subleading term proportional to the quantum depletion of the condensateu 

Let us further elaborate on the significance of this result. The depletion 
factor is intimately related to the fact that the Fock vacuum of the funda- 
mental bosons and the Fock vacuum of the phonons are inequivalent, a fact 
that is more profound that just an energy consideration. It is a statement 
about the full quantum state corresponding to the condensed phase, which is 
the regime in which we can speak about a semiclassical gravity analogue. It 
is a quantity that encodes in a specific way the information that the system 
is in a BEC phase, and not just a generic many body state. 

This observation, taken alone, would be of little significance outside a 
condensed matter community. However, it assumes a different relevance 
if we move to the perspective of a quantum gravity model. There, part of 
the properties of the long range/low energy /continuum/semiclassical regime 
that will result in a gravitational theory will be due to the microphysics, 
on one hand, and to the particular regime or considered state that will 
encode the fact that the model is considered in the long range/low en- 
ergy/continuum/semiclassical limit. The knowledge of the microscopic dy- 
namics must then be supplemented by the specific information about the 

In a comment [2T] to P2], an alternative point of view has been proposed. There it 
has been argued that the difference between the zero-point energy and the actual value of 
the cosmological constant is smaller, the closer is the analogue system to reproduce GR. 
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considered state. In this perspective, a similar situation arises within loop 
quantum gravity models [22j [23], suggesting a BCS energy gap as an origin 
for the cosmological constant. When we say that the cosmological constant 
has to be computed in terms of the microscopic theory, we also imply that 
we have identified the state or the class of states that will correspond to a 
semiclassical state, and that the effective dynamics will be the outcome of 
the microscopic dynamics as well as of the state considered. 

Going back to our model, this result suggests a twofold interpretation. 
First, there is no a priori reason why the cosmological constant should be 
computed as the zero-point energy of the system, even when this energy 
is calculated correctly taking into account the corrections coming from the 
microphysics of the system. The computation of this constant must pass 
instead through the derivation of Einstein's equations emerging from the 
underlying microscopic system. Second, the value of A can be several orders 
of magnitude smaller than the total vacuum energy density, depending on 
the features of the fundamental structure from which the spacetime emerges. 
The fine-tuning problem would therefore become much less worrying. How- 
ever, no indication about the coincidence problem comes from the analysis 
of this analogy. 

Moreover, even in such a simple system, several different scales show up. 
For instance, the Planck scale computed with the fundamental dimensionful 
constants is very different from the scale at which Lorentz symmetry is 
broken (i.e. when the dispersion relation is no longer linear), which in this 
case coincides with the range of gravitational interaction. This is just a 
coincidence due to the extreme simplicity of the model, and we have argued 
that a spinor BEC might alleviate this problem. In general, these scales 
are functions of the fundamental scattering length and of the diluteness of 
the condensate. Most importantly, the energy density associated with the 
cosmological constant can be much smaller than the value calculated with a 
cut off at the Planck scale, being proportional to the diluteness parameter 
of the condensate. Of course, this result must not be strictly translated 
to the gravity side of the analogy. It shows nevertheless how the problem 
of defining scales for possible Lorentz violations in relation to the so called 
Planck units is very far from being trivial. Furthermore, a very wide range of 
options is left open, depending on the fundamental structure the spacetime 
is emerging out. In particular, since the value of the cosmological constant 
appears to be strongly dependent on the granular structure of the spacetime, 
the comparison with the observed value of this constant might represent an 
important test for the validity of any theory of quantum gravity. 

Actually, we can further broaden the scope of the discussion. Our result 
strongly supports a picture where gravity is a collective phenomenon in a 
pregeometric theory (for a review of different ideas and references, see [24j ) . 
In fact, the cosmological constant puzzle is more easily addressed in those 
scenarios. From an emergent gravity approach, the low energy effective 
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action (and its renormalization group flow) for gravity and matter fields is 
obviously computed within a framework that has nothing to do with QFT 
in curved spacetime. 

In these scenarios, if we interpreted the cosmological constant as a cou- 
pling constant controlling some self-interaction of the gravitational field 
(i.e. putting it on the LHS of Einstein equations), rather than as a vac- 
uum energy (on the RHS), it would straightforwardly follow that the ex- 
planation of its value (and of its properties under renormalization) would 
naturally sit outside the domain of semiclassical gravity. From this point 
of view, the cosmological constant works as a phenomenological parameter 
controlling the constitutive relation that controls the self-interaction of the 
gravitational field as it is induced from (and at the same time it summarizes) 
the underlying dynamics. As such, its value cannot be computed within the 
semiclassical gravity regime, but has to be computed on the basis of the 
underlying microscopic dynamics. 

In this respect, it is conceivable that the very notion of cosmological con- 
stant as a form of energy intrinsic to the vacuum is ultimately misleading. 
This is the point of view of many quantum gravity models. An interest- 
ing case is represented by group field theories, a generalization to higher 
dimensions of matrix models for two dimensional QG [25]. To date, little is 
known about the macroscopic regime of them, conjectured to be some form 
of gravitational theory even though some preliminary steps have been re- 
cently done [26] . In such a framework, it is transparent that the origin of the 
gravitational coupling constants has nothing to do with ideas like "vacuum 
energy" or statements like "energy gravitates" , because energy itself is an 
emergent concept. Rather, the value of A is determined by the microphysics, 
and, most importantly, by the procedure to approach the continuum semi- 
classical limit. For instance, in [57] it is shown that the effective dynamics 
will be controlled essentially by the critical behavior of the model, as ex- 
pected, up to the point in which the gravitational coupling constants for an 
effective field theory will be functions of the critical exponents associated 
to the phase transition defining the continuum limit, i.e. to the state of the 
system that defines the continuum limit. Clearly these are concepts totally 
disconnected from notions like energy, fields, metric, etc. Most importantly, 
these critical exponents cannot be computed with field theoretic arguments 
in semiclassical gravity. The similarity to the case we are investigating here 
is manifest. 

Lacking a full fledged derivation of semiclassical continuum gravity from 
a quantum gravity model, toy models are playing a crucial role. While a 
detailed calculation cannot be avoided forever, analogue models can and will 
represent key assets, as sources of inspiration and techniques to address basic 
problems in quantum gravity. The simple model discussed here, involving 
a BEC in which an effective analogue dynamics for a gravitational field can 
be defined, shows how powerful this idea can be in elucidating the nature of 



22 



certain puzzles, as well as concretely suggesting a possible way to approach 
the solution in quantum gravity settings. In this respect, the reasoning 
of this chapter sheds a totally different light on the cosmological constant 
problem, turning it from a failure of EFT to a possible window on the process 
with which spacetime arises as an effective description. 
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